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ABSTRACT 

High resolution X-ray spectroscopy with the diffraction gratings of Chandra and XMM- Newton offers 
new chances to study a large variety of stellar coronal phenomena. A popular X-ray calibration 
target is Capella, which has been observed with all gratings with significant exposure times. We 
gathered together all available data of the HETGS (155 ks), LETGS (219 ks), and RGS (53 ks) for 
comparative analysis focusing on the Ne IX triplet at around 13.5 A, a region that is severely blended 
by strong iron lines. We identify 18 emission lines in this region of the HEG spectrum, including many 
from Fe XIX, and find good agreement with predictions from a theoretical model constructed using 
the Astrophysical Plasma Emission Code (APEC). The model uses an emission measure distribution 
derived from Fe XV to Fe XXIV lines. The success of the model is due in part to the inclusion of 
accurate wavelengths from laboratory measurements. While these 18 emission lines cannot be isolated 
in the LETGS or RGS spectra, their wavelengths and fluxes as measured with HEG are consistent 
with the lower resolution spectra. 

In the Capella model for HEG, the weak intercombination line of Ne IX is significantly blended by iron 
lines, which contribute about half the flux. After accounting for blending in the He-like diagnostic 
lines, we find the density to be consistent with the low density limit (ng < 2 x 10^" cm^"^); however, 
the electron temperature indicated by the Ne IX G-ratio is surprisingly low (~ 2 MK) compared with 
the peak of the emission measure distribution 6 MK). Models show that the Ne IX triplet is less 
blended in cooler plasmas and in plasmas with an enhanced neon-to-iron abundance ratio. 
Subject headings: atomic data - line: identification ~ stars: coronae - stars: individual (Capella) — 
stars: late-type - X-rays: stars 



1. INTRODUCTION 

Investigation of stellar coronae in the X-ray wavelength 
band has until very recently been restricted to instru- 
ments with intrinsically low spectral resolution, such as 
the proportional counters and CCDs on satellites such as 
Einstein, ROSAT, ASCA, and BeppoSAX.^ The hmita- 
tions of low resolution spectra have restricted X-ray stud- 
ies of stellar coronae to measurements of luminosities, 
plasma temperatures and estimates of elemental abun- 
dances. The spatial resolution of X-ray emitting plasma 
that is routinely available for studies of the solar corona 
with satellites such as SOHO and TRACE is still today 
a dream of X-ray astronomers. 

Despite limited spectral information, substantial 
progress has been made in understanding the gross char- 
acteristics of stellar coronae throughout the HR diagram. 
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^ Exceptions are the Einstein spectra of Cap ella with much bet- 
ter resolution, obtained by Mcwc ct al] 119821) with the objective 
grating spectrometer (OGS), covering the range 5 to 30 A with a 
resolution <1 A, and of Vc dder Canizares, 1.198,3.) with the focal- 
plane crystal spectrometer. 



For example, iSchriiver. Mewe. fc Walter! ^^1^84^ showed 
from Einstein IPC observations of a sample of 34 late- 
type stars that coronal temperature was directly corre- 
lated with X-ray luminosity and stellar rotation rate. 
Similar result s were obtained for a larger sample by 
iSchmitt et al.l l)199(](l . and later based on the so-called 
hardness ratio obtained from ROSAT PSPC observations 
by Schmitt (1997). 

These results raised a question as to the nature of 
the high tempe rature plasma in active stellar coronae. 
iVaiana & R.osn er ( 1978) had pointed out that the Sun 
completely covered with active regions would have an X- 
ray luminosity of ~ 2 x lO^^ergs"^. However, the most 
active solar-like stars can have X-ray luminosities up to 
two orders of magnitude higher than this and so cannot 
simply be scaled-up versions of the solar corona. The 
hot plasma on the most active stars must be structured 
differently from that in typical solar active regions. The 
radiative loss of a hot, optically-thin, collision-dominated 
plasma is essentially proportional to the volume emis- 
sion measure, defined as the product of electron density 
squared and the emitting volume, . Increasing ei- 
ther the emitting volumes or plasma density results in 
an increase in X-ray luminosity. 
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A key plasma parameter for inferring the size of X-ray 
emitting regions is therefore the plasma density. Direct 
spectroscopic information on plasma densities at coro- 
nal temperatures on stars other than the Sun first be- 
came possible with the advent of "high" resolution spec- 
tra (A/AA ~ 200) obtained by EUVE that were capable 
of separating individual spectral lines. Even with this 
resolution, the available diagnostics have often tended to 
be less than definitive, owing to the poor signal-to-noise 
ratio of observed spectra or to blended lines. This has 
been especially the case for the active stars. Studies of 
density-sensitive lines of Fe XIX to Fe XXII in EUVE 
spectra of RS CVn stars revealed tempting evidence for 
high densities of rie ^ 10^^ to 10^'^ cm~'^ at coronal tem - 
peratures near 10^ K l)Drakelll99d iDupree et al.lll993j) . 
Such high densities suggest that emitting structures are 
compact: static loop models suc h as those described by 
iRosner. Tucker, fc Vaianal l)1978f) would have heights of 
~ 1000 km with confining field strengths of up to IkG 
and surface filling factors of 1 t o 10%. In the case of 
the evolved active binary Capella. iBrickhou sc (1996) ob- 
tained tie ~ 10^2 cm-3 from Fe XIX to Fe XXII, but 
Tie ~ lO^cm"'^ based on Fe XII to Fe XIV, suggesting 
that the cool and hot plasma are from distinctly differ- 
ent structures. In contrast to the high densities reported 
by Dup rcc et al. (1993) only upper limits were found by 
iMewe et al.l l|20(yi^) (rie < 2 - 5 x 10^^ cm'^) using the 
same line ratios of Fe XX to Fe XXII but with the higher 
spectral resolution of Chandra LETGS spectra. They 
point out that the Fe XIX to Fe XXII line ratios are only 
sensitive above 10^^ cm"'^, such that no tracer for low 
densities for the hotter plasma component is available. 

The high spectral resolution of Chandra and XMM- 
Newton has significantly changed the situation regard- 
ing plasma diagnostics for stellar coronae. High spec- 
tral resolution coupled with large effective area allows 
application of line-based diagnostic techniques at X-ray 
wavelengths. Line ratios of the strong H-like Lyman se- 
ries, which are sensitive to the Boltzmann factor, and the 
He-like triplets, which exploit the competition between 
collisional excitation and recombination-driven cascades, 
can be useful temperature diagnostics. Simply examin- 
ing the presence or absence of principle lines of differ- 
ent ionic states of different elements gives an indication 
of the emitting plasma temperatures. The He-like sys- 
tems also provide density diagnostics based on the low- 
lying metastable level ls2s'^Si above the Is^ ^Sq ground 
state. The He-like oxygen diagnostic has been used to 
support the previous EUVE result for Capella, show- 
ing that the lower temperature plasma MK) also 
has lower densit y ([A udard et al. 2001; Brickhouse 200^ 
Brinkman et al.ll200(]l: iCanizares et al.ll200a INess et all 
200lir 

We compare the capabilities of this current generation 
of high spectral resolution X-ray instruments with a de- 
tailed study of the Ne IX triplet as a density and tem- 
perature diagnostic. We focus our study on the Capella 
binary system (HD 34029; a Aurigae; G8 III + Gl HI). 
Capella has been extensively studied from X-ray to radio 
wavelengths, and is the brightest, steady coronal source 
in the X-ray sky. As such, it has been a key calibration 
target for the Chandra high-energy and low-energy trans- 
mission grating spectrometers (HETGS and LETGS, re- 
spectively), as well as the XMM-Newton reflection grat- 



ing spectrometers (RGS), and has b een observed on sev- 
eral occasions by both satellites llAudard et alJ l20flll: 
iBrinkman et alJl2000t ICanizares et al.ll2000|) . The emis- 
sion measure distribution of Ca pella shows a steep but 
narrow enh ancement at 6 MK fBric khoiise et all 1200(1 
iDupree elTal. 1993) , making it ideal for studying the 
blending of Ne IX with high temperature lines. 

The purpose of this study is threefold: (i) assess the 
accuracy and reliability of state-of-the-art plasma radia- 
tive loss models for describing the emission spectrum of 
Capella in the region of the He-like complex of Ne; (ii) 
determine how well these models might describe the spec- 
tra of stars both more and less active than Capella and 
with different coronal temperatures; and (iii) gain further 
insight into the plasma density in the Capella coronae in 
the key temperature range 4 to 6MK. 

2. HE-LIKE IONS IN HIGH RESOLUTION STELLAR 
SPECTRA 

2.1. The Use of the He-like Triplet Diagnostics 

He-like ions are produced at temperatures ranging be- 
tween - 2MK for C V and N VI and 10 MK for Si XIII. 
The densities over which the diagnostics are sensitive in- 
crease with increasing element number Z, such that the 
lower Z ions C V, N VI, and O VII with temperatures 
of peak emissivity at 1.0, 1.6, and 2.0 MK, respectively, 
provide diagnostics for densities up to ~ 10^^ cm~'^. The 
higher temperature (> 6MK) ions Mg XI and Si XIII 
are sensitive at densities above ~ 10^^ cm"-^. Ne IX, 
typically formed at ^ 4MK, is sensitive to densities 
between ^ 10^^ and lO^'^cm"'^. In many cases mea- 
surements of Ne IX are the only chance to fill the gap 
between the cooler plasma and the hotter plasma. In 
jiddition, some s tars s h ow a neon overabundance (e.g., 
'Brinkman et a ll 120011: lOrake et all 120011: iGiidel et alJ 
l2001;:.Huenemoerder. Canizares fc SchulJ200lD. making 
Ne IX lines more easily detectable in those cases. How- 
ever, blending with Fe lines in the Ne IX triplet spectral 
region can compromise the diagnostic utility when ana- 
lyzing plasmas hot enough to produce them (especially 
Fe XIX at - 6MK). 

The t heory of the He-like triple ts was originally devel- 
oped bv lGabriel fc JordanI l)1969t) . The density diagnos- 
tic is often referred to as the i?-ratio, where R = f/i, 
denoting the forbidden line flux with / and the inter- 
combination line flux with i. The relation can be param- 
eterized as 

with low-density limit Rq and electron density 
rig. The critical de nsity Nr is the density at 
which R = \l2Rn. IPradhan fc Shulll llT98l and 
iBlumenthal. Drake, fc Tucked calculated theo- 

retical values of i?o and Nc and determined the Z- 
dependence of the density-sensitivity as shown in Fig- 
ure ^ Larger i?-ratios mean relatively weaker intercom- 
bination lines; since the i is intrinsically weak at low 
densities, it may be difficult to measure accurately. 

Strong ultraviolet radiation fields can also de-populate 
the metastable levels of the He-like triplets and change 
the i?-ratios. For the case of Capella and Ne IX, such 
fields can be neglected ( Ness et al....2001,.') . 

For an optically thin plasma, a temperature diagnostic 
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Fig. 1. — The low-density limits Ro and the critical densities 
A^c for the different He-like ions show systematic, anti-correlated 
trends with wavelength, due ultimately to the nuclear charge (from 
iniumcnthal. Drake, & Tucker 1972; Pradhan & Shull 1981). Large 
Ro (C V) implies that the i flux may be difficult to measure. Large 
Nc may be outside the interesting coronal range (Si XIII). Ne IX 
appears to have ideal values for coronal physics. 
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where r is the resonance line flux. In coUisional ionization 
equilibrium (CIE), the G-ratio decreases with tempera- 
ture primarily because the excitation of the triplet lev- 
els (by dielectronic recombination-driven cascades) de- 
crea ses faster than the collisional excitation to the ^Pi 
(see iSmith et al.l 1200 1|) . Assuming a priori that the 
plasma is in CIE, the G-ratio is a direct diagnostic of 
the temperature of Ne IX emission and can be compared 
to predictions based on an emission measure distribution. 

The He-like triplets have been used for measuring tem- 
peratures and densities in the solar corona. Both the 
O VII and Ne IX triplet s have been studied for quiescent 
emission and flares (e.g.. lActon et al.ll972HKeenan et al.l 
^ 87). In solar flares, the Ne IX lines appeare d blended. 
iMcKenziel (1985) and lDovle fc KeenanI l|198fiD suggested 
several candidate lines for these blends. 

Generally, stars with lower temperature coronae do not 
show significant blending around the Ne IX triplet. Pro- 
cyon (Ness ct al. 2002b) is a good example. In hotter 
coronae, the blending is more problematic. To date, the 
i?-ratios measured from MEG spectr a of Capella have 
not been used to determine dens ities l|Avres et al.ll200it 
iPhillips et aIll200lD . iNess et all (!2002a') have attempted 
to disentangle the blending for the LETGS/HRC-S spec- 
trum of Algol by assuming a priori a G-ratio of 0.8. 

2.2. Comparison among the New Instruments 

The new generation of X-ray telescopes of XMM- 
Newton and Chandra has opened new dimensions in sen- 
sitivity and resolution. Spectroscopic measurements can 
be carried out with both instruments using the intrinsic 
resolution of the CCDs (ACIS, EPIC) and the disper- 
sive instruments for higher resolution. Three gratings 
are providing data: the RGS on board XMM-Newton 
and the HETGS and LETGS on Chandra. The HETGS 



consists of two sets of gratings with different periods, the 
High Energy Grating (HEG) and the Medium Energy 
Grating (MEG), which intercept X-rays from the inner 
and outer mirror shells, respectively, and thus are used 
concurrently. The RGS and HETGS use the EPIC and 
ACIS-S CCD detectors, respectively, while the LETGS 
can use either ACIS-S or the micro-channel plate detector 
HRC-S. Figure |5] shows the effective areas and the wave- 
length ranges for the grating instruments. While RGS, 
HEG, and MEG operate in the spectral range 

^ 40 A, 

LETGS covers a much larger wavelength range from 2 to 
175 A. All the gratings have uniform resolving power for 
their entire wavelength range. 
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Fig. 2. — Effective areas for the Chandra and XMM-Newton 
gratings, with ranges of other observatories shown as dotted lines. 
At 13.6 A, the RGS has the largest area, while HEG is the least 
sensitive. The bottom panel shows detail in the Ne IX region. 



As can be seen from Figure |2 the large wavelength 
range covered by LETGS allows the extraction of total 
fluxes, luminosities, and hardness ratios corresponding to 
ROSAT (5 - 124 A) or Einstein (3 - 84 A). With LETGS 
all He-like triplets ranging from C V and N VI up to 
Si XIII can be measured simultaneously. LETGS also ob- 
tains An = Fe L-shell lines, including density-sensitive 
lines of Fe XIX to Fe XXII, previously obtained at five 
times lower spectral resolution with EUVE. 

Since LETGS cannot sort orders via detector intrinsic 
energy resolution, integrated fluxes and luminosities are 
still source- model-dependent. For a logT — 6.8 plasma. 
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the apparent energy flux obtained by dividing counts by 
the first order effective area is about 20% larger than the 
true first order flux, when 11 orders are accounted for. 
In this ideal isothermal model, only 80% of the counts 
are from first order. 

An advantage of the HETGS is its very high spec- 
tral resolution near 1 keV. In addition, order sorting for 
HEG, MEG, and the RGS is possible using the energy 
resolution of the CCD detectors, which is not possible 
for the LETGS/HRC-S configuration. This allows accu- 
rate determination of source-model-independent broad- 
band fiuxes and luminosities. 

The RGS covers a wavelength range similar to Einstein 
with especially large effective areas above 10 A {Fig.^. 
The spectral resolution of the RGS is comparable to that 
of LETGS in their region of overlap. 

The mirror point spread function primarily defines the 
dispersed line profile of the HEG, MEG, and LEG. In 
the case of the RGS, significant scattering wings also 
arise from the reflection gratings themselves. Grating 
period variance, detector pixelization, and aspect recon- 
struction are additional factors. For Ne IX, the RGS 
offers the highest sensitivity, even in second dispersion 
order; MEG and LETGS have larger effective areas than 
HEG. The spectral resolution of the HEG, however, is 
the most important factor for a realistic assessment of 
the effects of blending. Here we explore the diagnostic 
utility for observations of Ne IX with each instrument. 

2.3. Spectral Models 

Atomic data are fundamental to interpretation of the 
observed spectra. We compare observations with models 
produced u sing the Astrophysical Plasma Emission Code 
ijSmith et al. 2001) Version 1.2.^ APEC incorporates col- 
lisional and radiative rate data appropriate for modeling 
optically thin plasmas under the conditions of coUisional 
ionization equilibrium. The code solves the level-to-level 
rate matrix to obtain level populations and produces line 
emissivities as functions of electron temperature and den- 
sity. APEC also calculates spectral continuum emission 
from bremsstrahlung, radiative recombination, and two- 
photon emission. 

We rely on APEC in clusion of th e H- and He-like 
atomic data discussed bv lSmith et al.l l)2001fl . the HUL- 
LAC calculations of D. Liedahl for the Fe L-shell ions 
with additional R-matrix calcul ations available from 
CHIANTI V2.0 l|Landi et alJl999D . and isosequence scal- 
ing for Ni L-shell ions. Our models assume the ionization 
balance of iMazzo^l^ a et alJ lll998f) a n d the solar abun- 
dance model of lAnder^rGreve ssc (1989). Reference 
wavelengths are from quant um elec trodynam ic calcula- 
tions for H- and He-like ions l|Drakeiri9 88: Eri cssonll977fl 
and from laboratory measurements for emission lines of 
Fe L-shefl {Brown et al. 1998, 2002) and Ni L-sheU ions 
ijShira i et al. 2000). Additional L-shell wavelengths are 
derived from HULL AC energy levels. 

Despite significant improvements to spectral model- 
ing over the past several years, the theoretical atomic 
data in APEC and other plasma models remain largely 
untested over the broad range of applicable densities 
and temperatures. The deep observations of three late- 
type stars (Capella, Procyon, and HR 1099), used for 

^ Available at |http : / / cxc . harvard .edu / atomdb | 



in-flight calibration measurements of the dispersion rela- 
tion and line spread functions of the gratings, are also 
useful for determining the extent of agreement between 
models and observations and for highlighting issues that 
might require additional atomic physics work. The ac- 
curacy of each rate as well as the completeness of line 
lists is important to the correct interpretation of spec- 
tral diagnostic s, especially at lower spectral resolution 
(e.g., Brickhou se et alJ200(lD . This work is part of a com- 
prehensive effort known as the "Emission Line Project" 
to benchmark the atomic data in plasma spectral mod- 
els (Brickhouse & Drake! 120001) . Work to date on the 
Capella HETGS spectra has focused on the identification 
of strong lines for which HEG a nd MEG give good agree- 
ment (e.g., Avrcs ct al. 200^ iBehar. Cottam. fc KahnI 
120011 ICani zarcs et al. 2000J, while we aim to identify 
weak lines and test a comprehensive model. 

Atomic data for the He- like diagnostic lines themselves 
are exceptionally good — many rates have been bench- 
marked in controlled laboratory experiments. Toka- 
mak data confirm the calculation of the Ne IX R- 
ratio (Co ffev et al.lll99^ . The transition probability 
for the forbidden line, which determines the density- 
sensitivity of the Ne IX i?-ratio, has been measured 
to 1% accuracy on an electron beam ion trap (EBIT) 
fWargelin. Bcicrsdorfcr. & Kahn 1993); it has also been 
confirmed by B ragg crystal spe ctrometer measurements 
on the EBIT (IWargelinI 11993) . Smith et al. (2001) 
showed that the largest systematic uncertainties in the 
theory for O VII come from the limited number of energy 
levels used in calculating the cascades following dielec- 
tronic recombination. 

APEC models for Ne IX lines are in good agreement 
with the laboratory measurements. Simplifying assump- 
tions often found in the literature, such as the lack of 
temperature-sensitivity in the i?-ratio, are no longer nec- 
essary since the APEC models are calculated over dense 
temperature and density grids. 

For weak diagnostic lines such as the Ne IX in- 
tercombination line, two challenging line measurement 
issues — determining the continuum level and assessing 
fine blending — are greatly aided by APEC. The APEC 
line list includes approximately 40 lines with reference 
wavelengths between 13.35 and 13.85 A. Furthermore, 
the APEC emissivity table contains ~ 1000 lines in that 
region, down to 4 orders of magnitude weaker than the 
strong resonance line of Ne IX. This list includes Ne IX 
dielectronic recombination satellite lines and Fe and Ni 
L-shell lines with principal quantum number n < 5. 
APEC is sufficiently complete over the HETGS bandpass 
that it can be used to select line-free spectral regions for 
continuum fitting. 

3. OBSERVATIONS AND DATA REDUCTION 

Our approach is to use the six observations of HEG cal- 
ibration data of Capella to obtain a very long total expo- 
sure time (154.685 ks) and benchmark the APEC mod- 
els in the wavelength region around 13.6 A. The Capella 
observations discussed in this paper are summarized in 
Table n] For each of the three grating instruments, the 
data from several observations have been combined to 
produce high signal-to-noise spectra. Standard pipeline 
pr ocessing for HETGS, L E TGS, and RGS are des cribed 
in iCanizares et all (j2000(l . iBrinkman et alJ l|200H) . and 
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lAudard et al.l l)200H) , respectively. 

The HETGS/ ACIS-S data were obtained from the 
Chandra archive and reprocessed with CIAO software 
Version 2.2^ and calibration database CALDB Version 
2.10. The standard technique for separating the differ- 
ent grating orders employs a variable-width CCD pulse- 
height filter,"' that follows the pulse-height versus wave- 
length relation for each spectral order. Since some obser- 
vations were made during times of uncertain or changing 
detector response (primarily CCD temperature changes 
that affected the gain), the standard filter regions did 
not trace the actual distributions of events in the pulse- 
height-wavelength plane very accurately. To mitigate 
this, we used the constant fractional energy width op- 
tion available in the CIAO event-resolution program 
(tg_resolve_events) with a fractional width of 0.3. The 
constant width region accepts more inter-order back- 
ground, particularly at short wavelengths, but this is 
negligible for our purposes because of the very low back- 
ground level. We have also bypassed the step of po- 
sition randomization of events detected within a given 
pixel, which is part of the standard image processing of 
ACIS-S. Otherwise, we use the standard CIAO pipeline 
extraction. We use only the first dispersion orders. 
The effective areas for positive and negative first orders 
have been computed for each dataset using observation- 
specific data to account for aspect and bad pixels. The 
HEG effective areas for the individual observations in the 
wavelength range around 13 A are shown in Figure 
along with the exposure-weighted average effective area 
to be used for analysis of the combined spectra. 

The LET GS observa t ions a nd data reduction are de- 
scribed in iNess et a\} 1)200 ll^ . They consist of nine 
datasets added with a total exposure time of 218.54ks. 
Effective areas are from 2001 February (see Pease et al. 
2000 for a preliminary description). The RGS obser- 
vations were taken on 2000 March 25, with an expo- 
sure time of 5 2.92 k s. Data processing is described in 
lAudard et all l|2003|) . Effective areas have been calcu- 
lated for the observations with SAS5.3. 

We derive X-ray luminosities given in Table of Lx = 
1.3 X lO^o, 2.8 X lO^o, and 1.7 x lO^o 
distance of 12.94 ± 0.15 pc 



IjPerrvmi 

ferences are consistent with the different passbands and 
responses of the instruments. Table^also shows that the 
luminosities for HETGS and LETGS are fairly constant 
with time (for HETGS over 1.5 years), supporting our 
comparison of spectra taken with different instruments 
at different times. 
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4. DATA ANALYSIS 

4.1. Line Flux Measurements 

Table 121 lists 18 strong emission line features measured 
in the plus and minus first order HEG spectra between 
13.35 and 13.85 A, including several features not previ- 
ously reported for Capella. Line counts ar e measured 
with t he CORA line fitting tool deve loped bv INess et al.| 
( 20W ) and described and refined by INess &: W ichmann| 
1 2002(1 . The raw counts obtained with CORA represent 
the fitted number of expected counts for a given con- 

^ Available at http://cxc.harvard.edu/ciao 

* Encoded in the calibration database as the Order Sorting and 
Integrated Probability, or OSIP file. 
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Fig. 3. — Effective areas for the different observations (labeled 
with observation identification) with the HEG plus first order (top) 
and minus first order (bottom) . The average is indicated by a thick 
line and is obtained by weighting with the exposure times. The 
large dip on the short wavelength end of the minus first order is 
the chip gap, while smaller dips are caused by the removal of bad 
pixels. The widths and shapes of these dips are determined by the 
aspect dither and the chip geometry. 



tinuum (plus background), requiring Poissonian statis- 
tics to be conserved. Measurement errors are given as 
1 (7 errors and include statistical errors and correlated 
errors in cases of line blends, but include no system- 
atic errors from, e.g., the continuum placement. Line 
widths are fixed rather than fit, since line broadening is 
expected to be predominantly instrumental and line pro- 
files are represented by analytical profile functions (Gaus- 
sians for HEG and MEG, Lorentzians for RGS, and for 
the LETGS modified Lorentzians F(A ) = a/(l + ^^)'^ 
with (3 = 2.5: lKashvaD fc Drakd l|2002D V The continuum 
is chosen as an initial guess to be a constant with level 
determined by the HEG region at 13.6 A, which contains 
no apparent lines in the data as well as no significant 
emission lines in the APEC database. (The placement of 
the continuum is discussed further in ^.31 ) 
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Table 1. List of Observations 
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9nnn rr^ (t^ti fi-9^^-f^'? 

ZUUU-UO-UO -L -LO.ZO.OO - 


9000 0*^ OATni-l^i 


03 


30.10 


30.19 


01010 


29.541 


10773 


42907 


9nni 09 11T19-99-AQ 


9001 09 11T91-90 


17 


30.06 


30.12 


01099 


14.565 


5934 


23634 


1 QQQ nS 9S'Tn7-P^'^-1 7 


1 QQQ OS 9SX19-1fi 

iyyy-uo-zo -l iz. lo 


37 


30.12 


30.20 


01103 


40.479 


18109 


73303 


1999-09-24T06:09:21 - 


1999-09-24T18:22 


56 


30.16 


30.24 


01235 


14.571 


5916 


24177 


1999-08-28T12:16:37 - 


1999-08-28T16:39 


57 


30.12 


30.21 


01318 


26.701 


11133 


46338 


1999_09-25T13:26;39 - 


1999-09-25T21:51 


17 


30.14 


30.22 


HETOS Sum 


154.685 


63579 


257908 








30.12 


30.20 


Chandra LETCS 




LEC 










LEC^ 




01167 


15.36 


41215 




1999-09-09T13:10:06 - 


1999-09-09T17;26 


08 


30.45 




01244 


12.37 


32981 




1999-09-09T17:42:27 - 


1999-09-09T21;08 


36 


30.45 




01246 


15.00 


41701 




1999-09-10T03:06:06 - 


1999-09-10T07;16 


08 


30.47 




01248 


85.36 


240047 




1999-11-09T13:42:24 - 


1999-11-10T13;25 


05 


30.43 




01420 


30.30 


83533 




1999-10-29T22:49:29 - 


1999-10-30T07;14 


27 


30.44 




62410 


11.33 


31491 




1999-09-09T23:43:57 - 


1999-09-10T02;52 


48 


30.48 




62422 


11.68 


31202 




1999-09-12T18:26:42 - 


1999-09-12T21;41 


20 


30.46 




62423 


14.80 


39328 




1999-09-12T23:37:44 - 


1999-09-13T03;44 


28 


30.45 




62435 


22.34 


58201 




1999-09-06T00:35:40 - 


1999-09-06T06;48 


01 


30.42 




LETCS Sum 


218.54 


599699 










30.41 




XMM-Newton 




RCSl 


RCS2 








RGSl 


RCS2 


0121920101<i 


52.92 


169423 


186267 


2000-03-25Tll:36:59 - 


2000-03-26T02;53 


49 


30.22 


30.30 


0121920101= 


52.92 


45115 


44524 


2000-03-25Tll:36;59 - 


2000-03-26T02;53 


49 


30.29 


30.15 



''Counts are background subtracted values. The two dispersion directions are co-added. 
^Lx was determined within the wavelength range 3 to 20 A. 
'^Higher order photons are included without correction. 

First order. 
■^Second order. 



Table 2. Best-fit Line Fluxes for Separate Dispersion Orders of HEC 



A 






Amplitude'^ 


Flux'' 




Aoff 




A 






Amplitude'' 


Flux'' 


Acff 


(A) 






(cts) 


(ph cm ^ ks ^) 


{cm2) 


(A) 






(cts) 


(ph cm^^ ks^-') 


(cm^) 


-1-1^* order 












-l'^' order 










13.449 ± 





002 


315.9 ± 18.4 


0.511 ± 


030 


4.00 


13.449 


± 





003 


300.6 ± 17.9 


0.474 ± 0.028 


4.10 


13.553 ± 





008 


101.7 ± 10.8 


0.169 ± 


018 


3.88 


13.555 







006 


89.9 ± 10.1 


0.148 ± 0.017 


3.92 


13.700 ± 





004 


190.7 ± 14.4 


0.331 ± 


025 


3.72 


13.700 


± 





005 


182.8 ± 14.0 


0.312 ± 0.024 


3.78 


13.358 ± 





021 


31.1 ± 6.4 


0.050 ± 


010 


4.07 


13.356 







020 


31.1 ± 6.4 


0.048 ± 0.010 


4.21 


13.381 ± 





007 


40.6 ± 7.2 


0.065 ± 


012 


4.04 


13.378 







007 


39.5 ± 7.0 


0.061 ± 0.011 


4.17 


13.404 ± 





007 


52.0 ± 8.0 


0.084 ± 


013 


3.99 


13.404 


± 





007 


60.2 ± 8.4 


0.097 ± 0.013 


4.02 


13.427 ± 





009 


64.0 ± 8.7 


0.103 ± 


014 


4.01 


13.426 


± 





007 


74.8 ± 9.3 


0.117 ± 0.015 


4.13 


13.469 ± 





004 


134.8 ± 12.5 


0.219 ± 


020 


3.98 


13.468 







004 


145.3 ± 12.9 


0.236 ± 0.021 


3.98 


13.507 ± 





002 


193.6 ± 15.0 


0.318 ± 


025 


3.93 


13.507 


± 





002 


204.0 ± 15.7 


0.330 ± 0.025 


4.00 


13.524 ± 





002 


316.0 ± 18.7 


0.522 ± 


031 


3.91 


13.523 


± 





001 


261.7 ± 17.4 


0.425 ± 0.028 


3.98 


13.650 ± 





008 


71.2 ± 9.1 


0.122 ± 


016 


3.78 


13.646 


± 





007 


70.8 ± 9.2 


0.119 ± 0.015 


3.85 


13.675 ± 





009 


55.8 ± 8.3 


0.096 ± 


014 


3.75 


13.673 







009 


42.1 ± 7.4 


0.073 ± 0.013 


3.75 


13.722 ± 





006 


54.7 ± 8.2 


0.096 ± 


014 


3.70 


13.722 


± 





008 


47.7 ± 7.7 


0.082 ± 0.013 


3.76 


13.740 ± 





008 


42.2 ± 7.3 


0.074 ± 


013 


3.68 


13.742 







008 


50.2 ± 7.8 


0.089 ± 0.014 


3.66 


13.779 ± 





004 


98.4 ± 10.9 


0.175 ± 


019 


3.64 


13.778 







007 


94.3 ± 10.5 


0.166 ± 0.019 


3.67 


13.797 ± 





003 


146.6 ± 12.9 


0.262 ± 


023 


3.62 


13.797 


± 





006 


117.1 ± 11.4 


0.205 ± 0.020 


3.69 


13.828 ± 





004 


185.8 ± 14.1 


0.335 ± 


025 


3.59 


13.826 


± 





003 


172.0 ± 13.8 


0.313 ± 0.025 


3.55 


13.847 ± 





007 


50.4 ± 8.0 


0.091 ± 


014 


3.57 


13.843 







006 


55.4 ± 8.5 


0.102 ± 0.016 


3.50 



''Measured line counts with la errors. The line widths are all 0.005 A and the assumed continuum level (including background) 
for each dispersion order is 350 cts A"'^. 
''The total exposure time was 154.7 ks. 



The wavelengths and line counts are Usted in Table [3 
along with fluxes determined using the total exposure 
time and effective areas given in the last column for each 
line. We find the spectra from both dispersion directions 
to agree reasonably well and use the plus and minus first 
order-summed spectrum for further analysis (Table ■ 
The summed spectrum is shown in Figure 01 along with 
the empirical model, i.e., continuum plus 18 lines, with 



best-fit centroids and fiuxes. 

To predict count spectra from the other instruments, 
we convolved the HEG empirical model with each instru- 
mental response. We allowed global wavelength shifts 
(dA) to account for different absolute dispersion calibra- 
tions, and a normalization correction (Seal) for differ- 
ences in effective area calibration. A minimization 
adjusts these two parameters to transform the model to 
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Table 3. Line Identification for HEG Measurements^ 



(A) 


(A) 


(A) 


Ion 


Transition 


Fluxota 
(pli cm~2 ks~^) 


(pli cm~2 ks"-"^) 


1 Q Q C^/l 
iv5.O04 


1 Q Q c; c 

io.ooo 


u.uuy 


re XVIII 


^3/2 




U.U4o±U.UU / 


U.Uz/0 


lO.O I 1 


lo.oou 


U.UU4 


Pc XX 


2s2jj* ■* Pr /, 




U.UDZ±U.UUo 


U.Uo I i 




13.390 




Fe XX 


V ^■P3/2 






0.0046 


13.401 


13.395 




Fe XVIII 


V ^^3/2 




0.091±0.009 


0.0459 




13.407 


0.004 


Fe XVIII 


2p5 ^P3/9 
^1/ of Z 


- 2s2p? ,^2p? /o3pq/o 




0.0238 




13.409 




Fe XX 


2s2p* ,o 

^j^y J 5/2 


^0 ^^-j^ /2 ^i-'ci/ J '-'"'5/ z 




0.0169 




13.418 




Fexx 


ZSZp ^3/2 






0.0108 


13.424 


13.423 


0.004 


Fexix 


n_4 3 p„ 


— Zp [ U )ou -TS 


O.llOiO.OlO 


0.0504 




13.431 




Fexxi 


2s2p3 






0.0080 




13.432 




Fexxi 


2o'' ^ Pi /o 
-'1/2 


V ^ y'J'^ ^3/2 




0.0042 


13.446 


13.447 


0.004 


Ne IX 


Is^ l5n 


- Is2p Pi 


0.492±0.021 


0.3978 


io.4DO 




U.UUo 


Fe XIX 


2p4 3P2 


- 2p3(2D)3d 


U.ZZoitU.UlO 


U. 1140 


±O.OU4t 




U.UUO 


Fe XIX 




- 2pi/22p^/23'^3/2 


U . ZO III U . U i 


U.ZUUo 




13.507 


0.005 


Fexxi 


2s2p3 


- 2s2p2^23s 




0.0579 


13.521 


13.518 


0.002 


Fexix 


2p'* 3P2 


- 2p3(2D)3d 3P>3 


0.471±0.021 


0.4425 


io.ooi 


lo.OOU 


U.UUO 


Ne IX 


1^2 ISo 


- Is2p 3P2 


U.l0oitU.Ui.Z 


U.UUzl 




13.551 


0.005 


Fe XIX 


2p4 3P2 


- 2pi/22p^/23d5/2 




0.0290 




lO.UOO 


U.UUO 


Ne IX 




- Is2p 3 Pi 




U.UOOJ- 




13.554 




Fe XIX 


2p4 3P2 


- 2Pi/22p3/23(i5/2 




0.0101 




13.558 




Fe XX 


2s2p4 4p^/2 


- 2s2p? , 2p3/23(i3/2 




0.0112 


13.645 


13.645 


0.004 


Fe XIX 


2p'* 3p2 


- 2p3(2D)3d 3p3 


U.llo±O.Ull 


0.0712 




13.648 




Fe XIX 


2p4 3 P2 


- 2p3(4S')3d 3 




0.0175 


13.671 


13.674*^ 






2p4 


- 2pi/22p^ ,„3d5/2 


n nQc;_Ln ni n 
U. UOiJ^U. u ±u 


0.0133 




13.675 




Fe XIX 


2p4 3Pi 


- 2pi/22p^ ,,3(^3/2 




0.0264 




13.683 




Fexix 


2p4 3P2 


- 2pi/22p3/23d3/2 




0.0208 


13.697 


13.699 


0.005 


Ne IX 


Is^ 1 Sn 


— 1 s2s 3 


0.322±0.017 


0.1845 


13.719 


13.732' 




Fexix 


2p'l 3pi 


- 2pi/22p^/23d5/2 


0.089±0.010 


0.0573 


13.738 


13.746 




Fexix 


2p'* 1^2 


- 2p3(2D)3d 1P3 


0.081±0.010 


0.0690 


13.775 


13.767 


0.005 


Fexx 


*S3/2 


- 2p2(3p)3s 4p5/2 


0.170±0.013 


0.0494 




13.779 


0.005 


Nixix 


2p6 


- 2p5(2p)3s iPi 




0.0676 


13.794 


13.795 


0.005 


Fexix 


2p4 3P2 


- 2pi/22p2^23(ig/2 


0.234±0.015 


0.1778 




13.795 


0.005 


Fexix 


2p4 ID2 


- 2p3(2p))3d 3p2 




0.0207 


13.824 


13.825 


0.002 


Fe XVII 


2p'3 iSo 


- 2s2p'^3p iPi 


0.325±0.018 


0.3058 




13.839 


0.005 


Fexix 


2p4 3P2 


- 2p2^22P3/23rf5/2 




0.0321 




13.839 


0.005 


Fexix 


2p4 ID2 


- 2p2/22p3/23d5/2 




0.0082 


13.843 


13.843 


0.006 


Fexx 


"^3/2 


- 2p2(3p)3s 4p3/2 


0.094±0.011 


0.0233 



"Only the strongest lines modeled are identified. In cases of line blends, weaker lines up to about 10% of the strong line flux are also 
listed. 

'^Measured HEG wavelengths have been corrected according to the scaling. 
•^Sources for reference wavelengths with errors are given in the text. 

"^Model fluxes are based on the emission measure distribution derived from the lines in Table |4] (shown in Fig. ItI. 

•^Brown et al. (1999) list a group of 5 lines ("021") with a central wavelength 13.676 ± 0.004, although APEC has not assigned this 
wavelength to any of the model lines. 
'APEC includes two relatively strong Ncix satellite lines which may contribute to a centroid shift. 



each instrument. Figure |31 shows the MEG and LETGS 
model and data, while Figure shows RGSl and RGS2. 
The instrumental line spread function width (a) was not 
a free parameter, since it is well determined by calibra- 
tion. The continuum (cord) was also not a free parameter 
but was set at the flux level determined from the HEG 
spectrum. We find excellent agreement between the ob- 
served counts and the transformed HEG model, within 
the systematic errors expected from calibration. 

4.2. Emission Measure Distribution from the HEG 
Spectrum 

In order to compare our measurements with models we 
construct a rough emission measure distribution using se- 



lected iron lines from ionization stages Fe XV to Fe XXIV 
(see Table EJ. We use the strongest line from each ion 
(except for Fe XVII) and use other lines to check the 
model for consistency. The EUVE line fluxes are given 
by Avrcs ct al. (2001) and are corrected fo r an interstel- 
lar co lumn density Nh = 1-8 x lO^^cm"^ l|Linskv et alJ 
Assuming collisional ionization equilibrium, we 
construct an emission measure distribution at low density 
(rie = 1.0 cm~^) using the APEC line emissivities e{T). 
Each emission line specifies an emission measure curve 
iTTci^ FluXo6s/e(r), which represents the emission mea- 
sure as if each temperature contributed the total emis- 
sion in the line. It is an upper limit to the final, self- 
consistent emission measure distribution determined by 
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Ness, Brickhouse, Drake, & Huenemoerder 



in 



O 
u 




13.40 13.50 13.60 13.70 13.80 
X (A) 

Fig. 4. — The HEG summed spectrum of Capella using plus and 
minus first order. In order to construct an empirical model, eigh- 
teen emission lines are used with Gaussian line widths cr = 0.005 A 
(equivalent to FWHM= 0.012 A). The continuum level (including 
weak lines and background, for combined plus and minus first or- 
ders) is chosen constant with 700ctsA~^. The binsize is 2.5mA 
and the exposure time is 154.7 ks. 



using the lower envelope of the ensemble of curves as the 
initial estimate in an iterative optimization to a solution 
which predicts all measured line fluxes. Figure shows 
the individual emission measure curves with the best-fit 
model. The figure also compares an earlier emission mea- 
sure distribution obt ained from simultaneou s EUVE and 
ASCA observations ijBrickhouse et alJl200Cl|) . 

4.3. Modeling the Continuum 

Up to this stage of the analysis, the continuum has 
been taken as a constant value empirically determined 
from the spectrum. Since the value of the continuum 
under the weak lines is cruc ial to t he proper derivation 
of the i?-ratio fseerBrickho usel2002j) . we discuss our con- 
tinuum modeling method in some detail. In this section 
we derive a formal result which turns out to be close to 
our initial estimate; in general, iteration to subtract the 
new continuum flux from lines might be necessary. 

In principle, background, weak lines and continuum 
can be treated separately in the analysis. The back- 
ground, which includes non-source events from the detec- 
tor, cosmic rays, as well as source photons redistributed 
by mirror scatter, detector aliasing (CCD readout streak- 
ing), and grating scatter, is more than a factor of 20 
lower than the source-model continuum we derive. The 
HETGS spectra are not background-subtracted and so 
the background is implicitly included in the continuum 
model. However, the HETGS background rejection is 
very high: background near 13 A is estimated^ to be 
less than 0.017 cts ks~^ arcsec"^, and the extrac- 
tion width is 5 arcsec. Hence, we only expect about 13 
cts per order in the summed HEG spectrum. Weak 
lines not directly measured (i.e., not listed in Table |3Jl 
are treated in ti5.2l as a source of systematic uncertainty 
to the diagnostic line ratios. 



00 mA 

cont = 2776.12 cts A"' 
dA = -6.1 mA 
Scal = 0. 96207 
X /N= 10.54 




14.0 



1000 



800 



600 



400 



200 



CT = 90.0 mA 
cont = 3884.32 cts A ' 
dX=15.8 mA 
Scal = 1. 14989 
xVn = 5.812 




13.2 



13.4 



13.6 



13.8 



14.0 



^ See the "P roposers' Observatory Guide," 
|http: //cxc.harvard .edu/proposer/POG 



Section 8.3, 



Fig. 5. — Best-fit model obtained from the HEG spectrum, 
scaled and overlaid on the measured spectra from MEG (top) and 
LETGS (bottom). The scaling parameters dA, Seal, a, and cont 
(described in the text) are listed in the upper right box. The good- 
ness of fit x'^/N is also given. For MEG the binsize is 5.0 mA and 
the exposure time is 154.7 ks. For LETGS the binsize is 10.0 mA 
and the exposure time is 218.5 ks. 



The emission measure distribution constructed in tl4.2l 
was used to predict both continuum and line fluxes, 
under the assumption of solar abundances. The tabu- 
lated APEC linelist includes only lines with emissivities 
greater than 1.0 x 10"^" ph cm'^ s~^; weaker lines are 
included as a "pseudo-continuum" and are added to the 
continuum emission spectrum from bremsstrahlung, ra- 
diative recombination, and two-photon emission. The 
resulting model continuum is flat between 13 and 14 A 
and is about 30% lower than the continuum level used 
in t l4.1l (an initial guess chosen to match at 13.6 A), a 
difference of only about 3 counts per emission line. 

This model continuum spectrum was then fit to a 
set of line-free spectral regions, defined as spectral 
bins for which the model line flux is less than 20% 
of the model continuum flux. Both HEG and MEG 
were flt jointly using the Sherpa package in CIAO 
with the Cash statistic appropriate fo r low-count bins 
ijFreeman. Doe, fc Siemiginowskall200lD . With the con- 
tinuum shape fixed by the emission measure distribution, 
only a single parameter for the continuum level, or "nor- 
malization," was allowed to vary. The resulting fit con- 
tinuum level was 6% higher than the model continuum 
level calculated assuming solar abundances. Models with 
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Table 4. HEG Iron Line Measurements Useful for Emission Measure 
Distribution 







Ion 


Transition 


Counts 




(A) 


(A) 








(ph cm~-^ ks"-"^) (cm^) 



284.15'' 284.16 Fexv Ss^ i^o - 3s3p ^P2 •■• 10.5*" 

335.41'' 335.41 Fexvi S1/2 - ^p'^ Pi/2 ■■■ 33. 9*= 

15.262=^ 15.261 Fexvii 2p^ ^Sq - 2p^{^P)M^Di 841.008 ± 29.68 1.128 ± 0.040 4.819 

15.015 15.014 Fexvii 2p6 i^q - 2p5(2p)3d ^Pi 2719.34 ± 52.71 3.354 ± 0.065 5.241 

17.098 17.096 Fexvii 2p'^ ^ So - 2p^l^ P)Zs ^ P2 523.964 ± 22.95 2.459 ± 0.108 1.377 

14.207=^ 14.208== Fexviii 2p'^ ^ Pi/2 - 2p'^{^ D)Zd'^ D^/2 1347.92 ± 37.92 1.331 ± 0.037 6.547 

14.208<= Fexviii 2p5 ap^^^ - 2pi/22p3/23ci5/2 

16.075 16.071 Fexviii 2p5 apg/j - 2p4(3p)3s "Pg/a 286.749 ± 17.28 0.868 ± 0.052 2.135 

14.670^^ 14.664 Fexix 2p'^ ^ P2 - 2p'^{^ D)?,s Ds 164.356 ± 14.01 0.190 ± 0.016 5.586 

16.110 16.110 Fexix 2s2p^ ^P2 - 2pi/22plj^2,pi/2 42.796 ± 7.22 0.136 ± 0.023 2.037 

12.847'' 12.846 Fexx 2p3 4^^^^ _ 2pi/22p3/23d3/2 372.794 ± 20.84 0.222 ± 0.012 10.875 

12.828 12.824 Fexx "ip^ ^Sz/2 - '^Pi/2'^P3/2'idi/2 248.643 ± 18.42 0.148 ± 0.011 10.882 

12.827 Fexx 2p3 45-3/2 - 2pi/22p3/23rf5/2 

12.822 Fexxi 2s2p3 ^Di - 2s2pi/22p3/23d5/2 

12.286" 12.284 Fexxi 2p2 3p(, _ 2p3d 482.550 ± 23.27 0.233 ± 0.011 13.402 

11.771" 11.770 Fexxii 2s'^2p'^P-^/2 - 2s'^'id'^D.i/2 177.567 ± 14.80 0.077 ± 0.006 14.901 



12.750 12.754 Fexxii 2s2p"' ^03/2 - 2s2pi/23s 90.500 ± 10.40 0.050 ± 0.006 11.610 

11.741" 11.736 Fexxiii 2s2p ^Pi - 2s3d 107.723 ± 12.10 0.047 ± 0.005 14.666 

11.176" 11.176 Fexxiv 2p'^P^/2 - ^d^D^/2 90.073 ± 11.32 0.032 ± 0.004 18.211 



"This line is used to construct the emission measure distribution. 
''These measurements are from EUVE (Ayres et al. 2001). 
'^Sum of both hnes in the blend are used in the model. 



stricter criteria for determining line- free regions give sim- 
ilar results. The continuum level we adopt between 13 
and 14 A thus seems robust. 

Since the emission measure distribution was deter- 
mined only from iron lines, its absolute value should be 
lowered by 6% to give agreement with the fit continuum 
level, requiring a 6% higher iron-to-hydrogen abundance 
ratio. One can draw the preliminary conclusion from 
this analysis that the iron-to-hydrogen abundance ratio 
is close to solar; determining accurate abundance ratios 
requires additional checks from other spectral features 
and systematic error analysis which is beyond the scope 
of this paper. 

5. RESULTS AND INTERPRETATION 

5.1. Identification of the Emission Lines 

Using our derived emission measure distribution, we 
predict the spectrum of Capella under the assumptions 
of solar abundances and low density. We convolve this 
model spectrum with the instrumental response, allow- 
ing only for the small, systematic, global shift in wave- 
lengths, and compare it to the observed HEG spectrum. 
The model fluxes are given in Table |31 along with the 
measured fluxes, and the measured spectrum with this 
model is shown in Figure |H1 While the model is cal- 
culated at Ue = 1.0 cm~'^, we have used APEC VI. 3 
density-dependent calculations to check the sensitivity of 
our model iron spectrum. For < 10* cm^'^, no effects 
larger than a few percent are found. A few weak lines 
in the model show significant density-sensitivity above 
10* cm"-^, but no observable lines show more than 20% 
change up to = 10^-^ cm^-^. 

Except for one nickel line, only iron and neon lines are 
present in our model of this spectral region. The most 



striking discrepancies are seen for the Fe XIX lines at 
13.462 A (flux underpredicted by ~ a factor of two), at 
I3.497A (flux of blend underestimated by ~ one third, 
wavelength discrepancy consistent with errors on labora- 
tory wavelengths), and around 13.72 to 13.75 A (wave- 
length discrepancies well within theoretical errors, no 
laboratory wavelengths available). Minor disagreement 
is seen around 13.4 A. The nickel line at 13.779 A is un- 
derpredicted by a factor of two, probably an abundance 
effect. Since the goal is to assess blending of the neon 
diagnostic lines, the overall level of agreement is a good 
indicator of our knowledge of the blending spectrum (es- 
pecially for Fe XIX). We note that: (1) all strong pre- 
dicted lines are accounted for in the observed spectrum; 
(2) the worst flux disagreement is a factor of two; (3) 
laboratory iron wavelengths agree with the observations 
within the errors, while theoretical iron wavelengths ap- 
pear to agree to within one resolution element of the 
HEG spectrometer. 

Since we have derived the emission measure distribu- 
tion from iron lines, we can allow the three neon diag- 
nostic line fluxes to vary with independent scaling fac- 
tors, A([r, I,/]), in order to obtain best-fit neon abun- 
dances and densities. These parameters are iterated 
with a minimization. After applying the scaling fac- 
tor A{r) — 1.03 for the Ne IX resonance line, we find 
Flnxmodei = 0.41 ph cm"2 ks~^ 

The intercombination line (with an adjusted flux 0.075 
ph cm^^ ks^^) is significantly weaker in the model than 
in the observed spectrum. This is of course explained 
by the contaminating lines of Fe XIX and Fe XX that 
contribute to the observed feature. TableOland FigurelHl 
show that the iron lines can contribute almost half of the 
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CT = 35.0 mA 
cont = 721 .530 cts A ' 
dA = -0.6 mA 
Scal = 1.221 16 
xVn=1 -893 
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a = 35.0 mA 
cont = 805.595 cts A ' 
dX = 5.00 mA 
Scol = 1. 05598 
xVn = 3.034 



o 
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13.2 13.4 13.6 13.8 14.0 

Fig. 6.— Same as Figure Hffor^lhe RGSl (top) and RGS2 (bot- 
tom). For RGSl the binsize is 5.442 mA and for RGS2 the binsize 
is 6.048 mA. The exposure time is 52.92 ks. 




Fig. 7. — Best-fit emission measure distribution derived from 
iron lines i n ionization stage s Fe X V to Fe XXIV (solid). The 
result from l^icfchouse et alj 120001) is also shown (thick dotted). 
Individual line emission measure curves (thin dotted) are shown 
for the lines in Table HI 



total flux measured at 13.55 A, such that care must be 
taken with deblending this line even for HEG spectra. 
The forbidden line has an adjusted flux 0.28 ph cm~^ 
ks-i. 



5.2. Treatment of Weak Line Contamination 

Lines which are too weak to be identified and mea- 
sured in the spectrum (i.e., lines not found in Table|31) are 
treated as sources of systematic error to the diagnostic 
line fluxes. The total weak line flux in the region between 
13.35 and 13.85 A is predicted to be 1.9 times the model 
continuum flux; however, unlike the flat continuum, this 
line emission is not randomly distributed. Furthermore, 
most of the stronger of the unidentified lines have refer- 
ence wavelengths with small estimated errors. Thus it 
is reasonable to compute the fluxes of weak lines within 
the diagnostic line profile to estimate the degree of con- 
tamination. The models indicate that the resonance line 
contamination is ~ 5%, while the forbidden and inter- 
combination lines each have ~ 3% contamination. (This 
is in addition to the contamination from the 50% blend- 
ing of the intercombination line, already discussed.) 

5.3. Density and Temperature Diagnostics with Ne IX 
for Capella 

The Ne IX line fluxes constrain electron densities and 
temperatures. When blends are neglected, the Ne-Fe 
blended measurement (Table OJ yields an i?-ratio of 
2.1 ±0.3, which indicates a density log(ne) = 11.6 ± 0.1. 
Instead, accounting for the blends in the intercombi- 
nation line based on the emission measure model (see 
Fig. ^ we derive the low-density limit: R = 3.95 ± 0.7 
resulting in log(ne) < 10.2. Figure El shows that the 
low-density limit is found at all relevant temperatures. 

The G-ratio is less affected by blending. The raw mea- 
surements result in G = 0.97 ± 0.07, while the deblended 
fluxes give G = 0.93 ± 0.09, which leads to T= 2.1 ± 
0.9 MK (Fig.ini). 

The temperature derived from the G-ratio is a factor of 
two lower than the temperature of maximum emissivity 
(4MK), but probably consistent given the uncertainties. 
However, it is inconsistent with expectations based on 
the emission measure distribution. Figure IT^ shows that 
~ 94% of the line emission in our model comes from tem- 
peratures above 4MK. This discrepancy is perhaps most 
easily interpreted as an abundance effect in an inhomo- 
geneous plasma, though could also be caused by other 
effects. Standard emission measure analysis is not valid 
unless the emitting plasma has uniform abundances; this 
could easily be incorrect for active binary systems. 

If we allow for this possibility in the case of Capella, we 
could conclude that the dominant 6 MK plasma predicted 
by the peak in the emission measure distribution must 
have a lower than solar neon-to-iron ratio. However, if 
this were the case, in order to preserve the approximately 
solar neon-to-iron abundance ratio indicated by our line- 
to-continuum analysis (the scaling factor Ar = 1.03 for 
the Ne IX resonance line), the 2MK plasma must have 
a much higher neon-to-iron ratio. 

Young et al. ( 2001) have used FUSE to show that 
the Fe XVIII emission formed at 6 M K is associated 
predo minantly with the G8 giant, while l.lohnson et alJ 
piO^ . using HST/STIS, have found Fe XXI emission 
at ~ 10 MK to be predominantly on the rapidly rotating 
Gl giant. These latter observations were concurrent with 
some of those included in this analysis. Earlier measure- 
ments with HST/GHRS had found each star contribut ing 
roughly half the Fe XXI emission l|Linskv et alJll998j) . 
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dX = 2.97 mA 
Scal= 1.00000 
A(Ne) = (1.05, 1.45, 1.54) 
X/N = 6.366 



I— 

O 
O 




13.40 13,50 13.60 13.70 13.80 13.90 

X (A) 



Fig. 8. — Measured spectrum of Capella obtained from HEG (histogram with error-bars), APEC individual predicted emission Hnes 
(dotted Unes) and their sum (heavy solid line). Predicted line fluxes were estimated using the emission measure distribution described in 
84.21 Binsize is 2.5 mA and exposure time is 154.7 ks. 
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Fig. 9. — Same as Figure 1^ in more detail, showing the Ne IX 
intercombination line with contamination. 



Our results suggest that a high neon-to-iron abundance 



ratio is associated with the Gl star. Such a result in 
the case of the more rapidly rotating star would not 
be unprecedented. As noted in i)2.1l several analyses 
of Chandra and XMM- Newton spectra of active stars in- 
dicate abundance ratios for neon-to- iron considerably in 
excess of the accepted solar value. iDrake et all ll2001h 
also noted similar results arose in earlier analyses of 
low resolution ASCA studies. If the Gl Hertzsprung 
gap giant has 3x solar neon-to-iron ratio in its corona 
and contributes only 25% at the emission measure dis- 
tribution peak (based on the FUSE Fe XVIII measure- 
ment), then the clump giant would need to have less than 
about 40% of the solar neon-to-iron abundance ratio. A 
self-consistent analysis is difficult without more stringent 
constraints on the emission measure distribution, espe- 
cially below 6 MK. However, though the solar FIP effect 
presents an existing observational framework for a plausi- 
ble neon-to-iron ratio significantly below the solar value, 
it seems coincidental that the abundance ratio of the G8 
clump giant would conspire to produce a global average 
neon-to-iron abundance ratio which agreed so accurately 
with that of the Sun. Furthermore, preliminary analyses 
of the clump giants 7 Tau and f3 Get do not indicate low 
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4 . i/i=3^5+0_67 

T=10MK 




6 8 10 12 14 



Fig. 10. — Determination of the plasma density from the R- 
ratio using / and i from the direct measurement (//i = 2.1 ± 0.2) 
and from the model accounting for line blends (//i = 3.9 it 0.7). 
Shaded areas represent Itr errors on the measured ratios. The 
upper limit for the density is logne < 10.2. APEC models for the 
density-sensitive curves are shown for T = 2.0, 4.0, 6.3, 8.0, and 
10 MK, with the curve denoting the 6.3 MK model corresponding 
to the peak of the emission measure distribution. The _R-ratio 
increases with temperature. 



1 .2 




0.4 - 



0.2 r 1 

6.4 6.6 6.8 7.0 7.2 

Log T (K) 

Fig. 11. — Determination of the plasma temperature from the 
G-ratio for the fluxes adjusted for line blending (dashed line). Er- 
rors on the measured ratio are represented by dotted lines. APEC 
models are given for 10^", 10^^, 10^^, and lO^^cm"'^, with the 
solid curve denoting the 10^'^ cm""^ (i.e., low-density) model. 



neon-to-iron ratios (Brickhouse & Dupree 2003, in prepa- 
ration; Drake et al. 2003, in preparation, respectively). 

While we cannot totally rule out line blending, the re- 
quired contamination would be far larger than expected. 
If the G-ratio were high because of a contaminated / 
hne, the contaminating line would be the fifth strongest 
iron hne in this region, and therefore, hard for us to have 
missed. The discrepancy might also arise from errors in 
the ionization balance. The shape of the emission mea- 
sure distribution, derived primarily from the ionization 
balance of iron, is subject to significant atomic data un- 



0.60 b i I I i i 




Log T (K) 

Fig. 12. — The fraction of Ne IX resonance line emission aris- 
ing from each temperature in the plasma (solid curve), derived by 
multiplying the emission measure distribution of Figure|7|with the 
fractional emissivity curve from APEC (dash-dotted). The inter- 
combination and forbidden lines have virtually the same depen- 
dences and are not shown. 



certainties l|Brickhouse. Ravmond. fc Smithlll995(l . One 
further possible effect concerns the breakdown of the fun- 
damental assumptions implicit in the coronal approxima- 
tion: that the plasma is optically thin and is in ioniza- 
tion equilibrium. A consistent treatment of the inhomo- 
geneities in the system (two stars with different coronal 
structures contributing to the emission) is needed in or- 
der to test these assumptions. We defer further study of 
this complicated issue to future work. 

For Capella we have used iron lines observed in HETG 
and EUVE to determine the emission measure distribu- 
tion. For observations of other stellar coronae, which 
may have less sampling of the iron ionization stages or 
lower signal-to-noise than we have for Capella, the con- 
struction of an emission measure distribution becomes 
more complicated. A different approach is suggested 
(e.g.. ISchmitt fc Nessl l2003). where the emission mea- 
sure distribution is constructed by use of the H-like Lya 
and He-like resonance lines, including a model distribu- 
tion of magnetic loops. 

5.4. Blending as a Function of Plasma Temperature 

As was noted in H2.ll spectra of stellar coronae which 
are dominated by plasma at significantly lower tempera- 
tures than that of Capella show considerably less blend- 
ing around the Ne IX features, because the populations 
of the blending ions (primarily Fe XIX) are small. It is 
of interest to examine this more quantitatively in order 
to determine the temperature regimes for which Ne IX 
can be easily used. 

Figure ^| shows the cooling functions for the neon 
and blending iron lines. It is apparent that, in plasmas 
with temperatures logT below ~ 6.8, the measurement 
of Ne IX will not be severely complicated by blending. 
In plasmas with higher temperatures, the blending will 
be very strong unless neon is significantly overabundant 
compared to iron. Large neon-to-iron, as well as neon-to- 
hydrogen, ratios se em comm on amon g active stars, suc h 
as in H R 1099 (.Brinkman et al. 200LlDrake et alJ200l|) . 
II Peg l)Huenemoerder. Canizares. fc Schuldl200H) . and 
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AR Lac l)Huenenioerder et al.ll200'^ . Generally, in all 
kinds of low-pressure plasmas with high densities but 
low temperatures the blending with Fe XIX is less se- 
vere, and thus disentangling the Ne IX lines should be 
fairly straightforward. 

6. CONCLUSIONS 

The Ne IX triplet is an important tool for estimat- 
ing plasma densities, not only for coronal plasmas as 
in Capella, but for plasmas in general in the density 
range between 10^° and lO^'^cm"'^. However, we have 
shown that the intercombination line of Ne IX is severely 
blended in plasma with temperatures logT ^6.8. This 
is the case in active coronae as well as in solar flares, 
such that deblending of high temperature lines is also 
important for solar flare diagnostics. 
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Fig. 13. — Total emissivities as functions of temperature for each 
ion emitting strong lines in the wavelength range between 13.35 and 
13.85 A. Calculations are APEC models for solar abundances and 
low density. Table |21 gives the wavelengths for the lines identified 
in this work. 

Through our detailed study of Capella spectra, we have 
found that APEC models are sufhciently accurate and 
complete that all the significant observed lines in the 
Ne IX region can be reasonably identified in the HEG 
spectra. Th e laboratory wavelengths of Fe XIX from 
IBrown et alJ lj2002|) provide a significant improvement 
to the accuracy over the wavelengths derived from the 
HULLAC energy levels. The model fluxes for these lines 
are also in good agreement with the observations. In the 
APEC models for Capella, the Ne IX intercombination 
line is significantly blended. Given the predicted flux of 
the blending lines, the f/i ratio is consistent with the 
low density hmit (rig < 2 x 10^° cm""^). 

Interestingly, we find that the temperature-sensitive 
G-ratio is inconsistent with the emission measure dis- 



tribution derived from iron in the sense that it indi- 
cates significantly cooler electron temperatures. Since 
the strong peak in the emission measure distribution at 
6 MK is likely pro duced predominantly by the G8 star 
ijYoung et alJl200llh . abundance differences between the 
two coronae at first seem the most likely explanation for 
this apparent inconsistency. This would require the neon- 
to-iron abundance ratio to be lower in the G8 corona than 
in the Gl corona. In contrast, the neon-to-iron ratio in 
the Gl corona would have to be higher than solar values 
such that the average ratio in observations of the Capella 
system as a whole appeared solar. 

If correct, this interpretation could have significant im- 
plications for coronal analyses in general, since many 
coronal sources are binary systems comprised of two coro- 
nae. Furthermore, if coronae on individual stars are com- 
positionally inhomogeneous (as is the case on the Sun), 
models will have to account for this. In future work, we 
will investigate more detailed models to determine what 
abundance and temperature differences are required to 
explain the Capella G-ratio and emission measure dis- 
tribution inconsistency. We will also continue the de- 
tailed assessment of blending and atomic data for these 
and other spectra, which may shed further light on this 
problem. 

In the case of Capella spectra obtained by LETGS and 
RGS, the spectral resolution is inadequate to derive the 
same results independently from this isolated spectral re- 
gion. We have not yet explored the possibility of using 
other stronger, isolated lines from Fe XIX to Fe XXI to 
specify the iron contribution to the blended Ne IX spec- 
tral region, and then perform a more constrained fit. The 
potential problem with this approach is that typical un- 
certainties in the APEC model line flux ratios (~ 20 to 
30%) may be too large to sufficiently constrain the blend- 
ing spectrum. With a single test case such as Capella, 
we cannot explore all the parameter space in the mod- 
els, and thus any conclusions concerning the treatment 
of blending cannot yet be generalized. Nevertheless, the 
good agreement between the APEC models and the HEG 
observations is encouraging. 
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